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Abstract 

The method of characteristics has played a very important role in mathematical 
physics. Preciously, it was used to solve the initial value problem for partial 
differential equations of first order. In this paper, we propose a fractional method of 
characteristics and use it to solve some fractional partial differential equations. 
Keywords: Modified Riemann-Liouville Derivative; Fractional Method of 
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1 Introduction 

In the past centuries, many methods of mathematical physics have been developed to 

solve the partial differential equations (PDEs) [1 - 2], among which the method of 

characteristics is an efficient technique for PDEs [3]. 

Fractional partial differential equations have attracted many researchers' interests. 
Then one of the questions may be naturally proposed: would it be possible to derive 
the exact solutions of partial differential equations (FPDEs) using a fractional method 
of characteristics? Recently, with the modified Riemann-Liouville derivative [4], G. 
Jumaire ever proposed a Lagrange characteristic method [5] which can solve the 
FPDEs 

a{x,t) \ ,J +b(x,t) \ ,J =c(x,t), 0<ot,p<l, (1) 

ox p at 

where a = p. However, in most cases, the two fractional order parameters a and 

P may not be equivalent. 

In this letter, we present a more general fractional method of characteristics and 

consider the case a ^ p. 



* Corresponding author, E-mail address: wuguocheng2002@yahoo.com.cn (G.C. Wu). 



2 Properties of Modified Riemann-Liouville derivative and Integral 
Through this paper, we adopt the fractional derivative in modified Riemann-Liouville 
sense. The modified Riemann-Liouville derivative has been successfully applied in 
successfully applied in the probability calculus [6], fractional Laplace problems [7], 
fractional variational calculus [8] and fractional variational iteration method [9]. 

Firstly, we introduce the definitions and properties of the fractional calculation. 

Assume / : R — > R : x — > f(x) denote a continuous (but not necessarily differentiable) 

function and let the partition h > of an interval [0, 1]. Through the fractional 
Riemann Liouville integral 

/;/(x) = -i- \\x-5r*f{&d$ t a > 0. (2) 

T{a) J0 

The modified Riemann-Liouville derivative is defined as 

d:/(x) = — ^— j- \*(x -#-(/(£) -/(0))rf£ (3) 

where xe[0,l], 0<a<l. One also can obtained Eq. (3) through consequence of a 
more basic definition, a local one, in terms of a fractional finite difference [10] 



r 
A a = fW- 1") / (x §=- ( 



n f+nt- (k (4) 



where FWf{x) = f(x + h). Then the fractional derivative is defined as the following 
limit 

/(a)=lim A^x) 

h->0 H 

Some properties for the proposed modified Riemann-Liouville derivative are 
listed 

(a) Fractional Leibniz product law [10] 

D ( f(uv) = u (a) v + uv ia \ (6) 

where u (a) = D ( "\u) . Much more generally, we can have 



D { x " a) (uv) = Y j Cy ka) v (( "- k)a) . (7) 

fc=0 

(b) Fractional Leibniz Formulation [10] 

Q r x D a x f{x)=f Ji-y (o),a<o< (8) 

Note that the properties (a) and (b) lead to the integration by parts: 

Xu (a \ = {uv)l b a -Xuv (a \ (9) 

(c) Fractional Jumarie-Taylor series [10] 

oo foak 

f( x + h) = Y J -—J (ak) (x). (10) 

t~i (ak)\ 

(d) Integration with respect to (dx) a 
Assume f(x) denote a continuous i? — > i? function. We use the following 
equlitiy for the integral w.r.t (dx) a [ 1 0] 

r(«) Jo r(«+i) Jo 



(e) Some useful formulas 



/(M0D Ca) =f-* (a) W; 

dx 

£)(<*) X P _ T(l + P) „P-a 



r(l + p-a) ' ( 12 ) 

I (dx) = x ; 
r(l + a)Jx = d a x. 



3 Fractional Method of Characteristics 

It is well known that the method of characteristics has played a very important role in 

mathematical physics. Preciously, the method of characteristics is used to solve the 

initial value problem for general first order. Consider the following first order 

equation, 

. . 8u(x,t) .. .du(x,t) , s „„. 

a(x,i) — —^- L + b(x,t) — ^- L + = c(x,f). (13) 

dx dt 



The goal of the method of characteristics is to change coordinates from (x, t) to 

a new coordinate system (x , s) in which the PDE becomes an ordinary differential 
equation (ODE) along certain curves in the x — t plane. The cures are called the 



characteristic curves which read 



du , 
as, 



c(x,t) 

dx 
a(x,t) 

dt 



ds, 
ds. 



b(x,t) 

With the modified Riemann-Liouville derivative, G. Jumaire ever gave a 
Lagrange characteristic method [5] which can solve some classes of fractional partial 
differential equations. In this paper, we present a more generalized fractional method 
of characteristics and use it to solve linear fractional partial equations. 

More generally, we extend this method to linear space-time fractional differential 
equations 



^ u(x, t ) d": u ( x,j) 



a(x, t] y ' 4- &(x,-0— — = J c( x, Qq. (S)< (14) 



Expand u as the fractional Jumarie-Taylor's series of multivariate functions [10], 



d u{x,t) p d a u{x,t) 

r(l + p)cbc p + T(\ + a)df 



du = £; v ;;; 7 p (dx) + Z. 1 „ ( dt T, < a, p < 1. (15) 



The total derivative here is more generalized. The function u here is (3"" order 

differentiable with respecter to x and a' h order differentiable to t , respectively. 

Similarly, note that the generalized characteristic curves can be presented by 

du , . 
— = c(x,t), 
ds 

(dx) . . 

— - — = a(x,t), (16) 

— ^^ = b(x,t). 

T(l + a)ds 



Eq. (16) can be simplified as Jumaire's Lagrange method of characteristics (See. 
A. 12 in Ref. [5]) if we assume a = |3 

(dxf _ (dtf _ T(l + a)du _ d a u 
a(x,t) b(x,t) c(x,t) c(x,t) 

Obviously, if a = 1 in Eq. (19), we can get the characteristic curve for Eq. (13) 

a(x, t) — ) -^- L + b(x, t) — ^-^ = c(x, t). 
dx dt 

4 Application of Fractional Method of Characteristics 

Example 1 . As the first example, we consider space-time fractional equations for 
the transport equation in porous media, 

d a u(x,t) d u(x,t) „ n _ . 

Assume Eq. (20) subjects to the initial value u{x, 0) = (p(x). 

The generalized characteristic curves satisfy 

du 

ds 

(dx) P =cT(l + P)ds, (19) 

(dt) a =T(l + a)ds. 



Then we can obtain 



x , n 



r(i+p) 

-f— = s + C 2 , (20) 

T(l + a) 

u = C 3 . 

where C x , C 2 and C 3 are integral constants. Eliminating the parameters, we find 



x ct a 
that the fractional curves = x n and h is a constant along the 

r(l + P) T(l + a) ° 
fractional curves. 

Then we can directly derive the exact solution of Eq. (18) has the following form 



x ct a 

u = f(— ), (21) 

T(l + p) T(l + a) 

X x ct 

where f( ) = cp(x). By setting u(x,t) = f( ), we directly 

T(l + p) T(l + p) T(l + a) 

have an exact solution for the initial-value problem here. 

Example. 2. As the second example, we investigate the more complicated 
equation, 

'- dau{XJ K X " d ^ t} =0, 0<a, P <l. (22) 

T(l + a) dt a T(l + p) 5x p 

We can have the generalized cure equations 



ds 




(dx) 


x p 


r(i+p)* 


r(i+p)' 


(dty 


t a 



(23) 

T(l + a)ds T(l + a) 
With the properties (12), direct calculation leads to 

= c 2 e s , (24) 



x p 

r(i+p) 

f 



T(l + a) 

u =c 3 . 

Then we can find Eq. (22) has the following solutions 

x p , f 



u=f( / ), (25) 

T(l + P) T(l + a) 



where the function / is arbitrary . 



x t a 

Assume X = / . With the properties (12), we note that 

T(l + P) T(l + a) 



O u(x,t) _ y(a) _ _ f X t 

8t a ~ x ' *r(i+p) (r(i+a)) 2 ' 



: f x X { , a) = ~fx ^-^ I „„' xx2 , (26) 

or ' - 

and 



dUiX ' t) =f x XT=f x l^^. (29) 



5x p '' ' x T(l + a) 

As a result, we can proof 



t a d a u(x,t) x p d u(x,t) 



r(i+a) a; a r(i+p) ax p 

~ ra+p) ra+a) ra+p) r(i+a) 

= 0. 

5. Conclusion 

The classical method of characteristics is an efficient technique for solving 
partial differential equations. In this paper, fractional method of characteristics is 
under consideration for some classes of fractional partial differential equations and 
two examples are illustrated to show its efficiency. Besides, the presented method 
provides a potential tool to solve fractional symmetry equations in Lie group method. 
We will discuss such work in future. 
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